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Abstract. We consider a class of branching-selection particle systems on R sim- 
ilar to the one considered by E. Brunet and B. Derrida in their 1997 paper "Shift 
in the velocity of a front due to a cutofF. Based on numerical simulations and 
heuristic arguments, Brunet and Derrida showed that, as the population size N 
of the particle system goes to infinity, the asymptotic velocity of the system con- 
verges to a limiting value at the unexpectedly slow rate {logN)~'^. In this paper, 
we give a rigorous mathematical proof of this fact, for the class of particle systems 
we consider. The proof makes use of ideas and results by R. Pemantle, and by 
N. Gantert, Y. Hu and Z. Shi, and relies on a comparison of the particle system 
with a family of A'' independent branching random walks killed below a linear 
space-time barrier. 



1. Introduction 

1.1. Brunet-Derrida behavior. In [6l [7], E. Brunet and B. Derrida studied, 
among other things, a discrete-time particle system on Z, in which a population 
of particles with fixed size undergoes repeated steps of branching and selection. 
As time goes to infinity, the population of particles, taken as a whole, moves 
ballistically, with an asymptotic speed depending on the population size N. One 
remarkable property of this system is the following: as goes to infinity, the asymp- 
totic speed of the population of particles converges to a limiting value, but at the 
unexpectedly slow rate of (log A^)~^, bringing to light an unusually large finite-size 
effect. This behavior was, on the one hand, observed by Brunet and Derrida on 
direct numerical simulations of the particle system (with large numbers of particles, 
up to = 10^^). On the other hand, Brunet and Derrida provided a justification for 
this behavior through the following argument. First, in the limit where A^ goes to 
infinity, the time-evolution of the distribution of particles in the branching-selection 
system is governed by a deterministic equation, which can be viewed as a discrete 
version of the well-known F-KPP equation 

(1) ^ = Au + u{l-u), 

where u = u{x,t), x G M, t > 0. To account for the fact that there is only a finite 
number A^ of particles in the system instead of an infinite one - whence a resolution 
equal to 1/A^ for representing distributions of mass -, one may introduce a cut-off 
value of 1/A'" in the equation, and expect that this modified equation still refiects at 
least some of the behavior of the original particle system. Whence the question of 
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studying, for large A^, an equation of the form: 

du 

(2) — = Au + u{l-u)l{u>l/N). 

In fact, Brunet and Derrida could provide heuristic arguments for this new problem, 
showing that, for large A^, the effect of the cut-off is to shift the speed of travelling 
wave solutions of Equation ^ from the speed of those of Equation ([1]), by an 
amount of order (log A^)^^. In turn, these arguments were supported by numerical 
simulations of (discrete versions of) Equation ([5]). This result concerning the F- 
KPP equation with cut-off has recently been given a rigorous mathematical proof, 
see [21E1III]. 

A related question (see [8]), is that of the behavior of the F-KPP equation with 
small noise, i.e. of the equation 

(3) — = Au + u{l-u) + y ^ ^ V , 

where is a standard space-time white-noise, and A^ is large. Rigorous results have 
recently been derived for this model too, see [9l[Ill[l5j, establishing that the speed 
of the random travelling wave solutions of Equation ([3]) is, for large A^, shifted from 
the speed of those of ([I]) by an amount of order (log A^)~^. 

We thus have (at least) three examples of what may be called Brunet-Derrida 
behavior, in three different and more or less loosely related frameworks (branching- 
selection particle systems, F-KPP-equation with cut-off, F-KPP equation with 
noise), two of which have already been established rigorously. 

1.2. Main result. The goal of this paper is to give a proof of Brunet-Derrida 
behavior for a class of branching-selection systems that is similar (but not exactly 
identical) to the one originally studied by Brunet and Derrida in [6l [7] . 

To be specific, we consider a discrete-time particle system with A" particles on R 
evolving through the repeated application of branching and selection steps defined 
as follows: 

• Branching: each of the A^ particles is replaced by two new particles, whose 
positions are shifted from that of the original particle by independently per- 
forming two random walk steps, according to a given distribution p; 

• Selection: only the A^ rightmost particles are kept among the 2A^ obtained 
at the branching step, to form the new population of A^ particles. 

Our assumptions on the random walk distribution p are listed below, and come 
from the need to apply the result of the paper [13] by N. Gantert, Y. Hu and Z. 
Shi on the survival probability of the branching random-walk killed below a linear 
space-time boundary, in the special case of deterministic binary branching. 

Introduce the logarithmic moment generating function of p defined by 

A(i) := log J exp{tx)dp{x). 

Here are the assumptions on p: 
(Al) The number a := sup{t > 0; A{-t) < +00} is > 0. 
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(A2) The number C := sup{t > 0; A{t) < +00} is > 0. 
(A3) There exists t* g]0,C[ such that t*A'{t*) - A{t*) = log 2. 
Under these assumptions, both numbers 

Xip) := 4t*A"it*), v{p) := A'(t*) 

are well-defined, and satisfy < x{p) < +00 and v{p) G R. Simple cases for 
which these assumptions hold are e.g. the Bernoulli case for a €]0, l/2[, where 
p = aSi + (1 — a)6o, the uniform case, where p is the uniform distribution on the 
interval [0, 1] , and the gaussian case, where p is the standard Gaussian distribution 
on M. 

In Section [3] below, it is proved that, after a large number of iterated branching- 
selection steps, the displacement of the whole population of N particles is ballistic, 
with deterministic asymptotic speed vn{p), and that, as N goes to infinity, viy{p) 
increases to a limit Voo{p), which turns out to be equal to the v{p) defined above, 
and is thus finite under our assumptions. The main result concerning the branching- 
selection particle system is the following theorem: 

Theorem 1. Assume that (A1)-(A2)-(A3) hold. Then, as N goes to infinity, 
(4) Vooip) - vn{p) ~ x{p){^ogN)-^. 

1.3. Credits. The proof of Theorem [1] given in this paper is based on a compari- 
son of the particle system with a family of N independent branching random walks 
killed below a linear space-time barrier, and makes use in a crucial way of ideas and 
results from the following two sources: the paper [16] by R. Pemantle on complexity 
bounds for algorithms seeking near optimal paths in branching random walks, and 
the paper [13] by Gantert, Hu and Shi on the survival probability of the branching 
random-walk killed below a linear space-time boundary. A detailed description of 
exactly which ideas and results are used and how is given in Sections U [5] and [6] 
below. Note that the existence of a link between the Brunet-Derrida behavior of 
a branching-selection particle system such as the one studied here, and the asymp- 
totics of the survival probability for branching random walks killed below a linear 
space-time barrier, was already suggested in the papers [lOi, il7J by B. Derrida and 
D. Simon, where Brunet-Derrida-like features were observed for a quasi-stationary 
regime of killed branching random walks; the present paper gives an explicit and 
rigorous version of such a relation. 

Finally, let us mention that a first version [3] of the present work was completed 
by one of the authors (J.B.) before the results in [13] became publicly available. 
In [4], only the (logA^)"^ order of magnitude of the difference Voo{p) — vn{p) in 
the Bernoulli case was established. The results in [13] then allowed us to prove 
Theorem [H which is both more precise and more general. 

1.4. Organization of the paper. The rest of the paper is organized as follows. 
In Section [21 we provide the precise notations and definitions that are needed in 
the sequel. Section [3] contains a discussion of various elementary properties of the 
model we consider. Section [4] collects the results from [L3l that are used in the sequel. 
Section [5] contains the proof of the lower bound part of Theorem [U while Section [6] 
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contains the proof of the upper bound part. Section [7] discusses the Bernoulh(a) 
case for a > 1/2, showing that the conclusion of Theorem [T] may fail to hold when 
Assumption (A3) is not met. Section [8] is an attempt to provide a self-contained 
explanation of the {logN)~'^ order of magnitude appearing in Theorem [H The 
arguments in this section are only discussed in an informal way. 

2. Notations and definitions 

2.1. Particle systems on M. It is convenient to represent finite populations of 
particles by finite counting measures on M. We use the notation C to represent the 
set of all finite counting measures on M. 

For v £ C, the total mass of u (i.e. the number of particles in the population 
it describes) is denoted by M[u). We denote by maxz/ and minz^ respectively the 
maximum and minimum of the (finite) support of u. We also define the diameter 
d{v) := maxz^ — mini^. Given /i, € C, we use the notation -< to denote the usual 
stochastic ordering: fi ^ v if and only if fj,([x, +oo[) < u{[x,+oo[) for all x G M. 
In particular, ^ ^ v implies that M{p,) < M(i/), and it is easily seen that, if 

= Yli=i^ ^x, and u = Ylii^i^ '^ith xi > ■ ■ ■ > xm{^) and yi > ■ ■ • > yM{u), 
^ ^ V \s equivalent to M{fj,) < M{u) and Xi < in for all i G p,M(/x)]]. 

For all > 1, let Cat denote the set of finite counting measures on M with 
total mass equal to A^. In the sequel, we use the notation {X^)n>o to denote 
a Markov chain on C^v whose transition probabilities are given by the branching- 
selection mechanism with A^ particles defined in Section \1.2\ and which starts at a 
deterministic value Xq S Cat. We assume that this Markov chain is defined on a 
reference probability space denoted by (0, J-", P). 

2.2. Branching random walks. In the sequel, we use the notation BRW to denote 
a generic branching random walk on a regular rooted binary tree, with value zero 
at the root, and i.i.d. displacements with common distribution p along each edge. 
More formally, BRW consists of a pair (T, <I>), where T is a regular rooted binary 
tree, and <I> is a random map, associating to each vertex u G T a random variable 
^{u) G M in such a way that $(root) = and that the collection {^{v) — $(ii))^^ 

is i.i.d. with common distribution p, where {u, v) runs over the set of pairs of vertices 
of T such that u is the father of v. We say that $(ti) is the value of the branching 
random walk at vertex u. The probability measure governing BRW is denoted by 
Q. 

Given m > 1, we say that a sequence uq, . . . , Um of vertices in T is a descending 
path if, for all i G [1, mj, Uj_i is the father of Ui. The set of vertices of T located at 
depth m is denoted by T(m). 

3. Elementary properties of the model 

As a first quite elementary property, note that, from Assumptions (Al) and (A2), 
E(maxA^) and E(min A^) are finite for all n > 0, for any choice of the (determin- 
istic) initial condition Xq G Cat. 
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3.1. Estimates on the diameter. 

Proposition 1. Let un' '■= r \og^ "l ~'~ ^' ^'^^ '^^^ N >1, all initial population Xq € 

Cat, and all n > un, d{X^) is stochastically dominated by un x (m^^ — mj^^), where 

wS^ and ni'j!^ are respectively the maximum and the minimum of a family of2Nu]\i 
i.i.d. random variables with common distribution p . 

Proof. Consider n > U]\f, y = maxX^„^, and let us study the evolution of the 

branching-selection system between times n—ujsi and n. Define ni'j!^ as the minimum 
of the 2Nun random walk steps performed by the system between times n — uj\f and 
n. 

Consider first the possibility that minX^ < y + {k — n -\- U]\j)m^^^ for all k € 
[n + l — utv, Since all the random walk steps that are performed during branching 

steps are > m^^\ this implies that all the particles descended by branching from a 
particle located at y at time n — ujy , are preserved by the successive selection steps 
performed from X^_^^ to X^ . Since there are at least 2'^'^ > N such particles at 
time n, this is a contradiction. As a consequence, we know that there must be an 
index G [n+l — li^r, n| such that min > y + {k — n+UN)m^^\ Again by the fact 
that random walk steps are > m^^\ 1 1— t- min X^ — {t — n-\-UN)m^^^ is non-decreasing 
on the interval [n -|- 1 — utv, "^1, so we deduce that minX^ > y + UNm!"^ . Now, let 

(2) 

m]y denote the maximum of the 2Nu]\f random walk steps that are performed at 
the branching steps between time n — u^ and time n. We see from the definition of 
y that maxX^ < y + UN'm^N • have also just seen that minX^ > y + unit^^n ^ 

T-N J- 

□ 



SO that d{X^) = maxX^ — minX^ < ui\]{m^'^ — rnf^^ 



The following corollary is then a rather straightforward consequence, in view of 
Assumptions (Al) and (A2). 

Corollary 1. For all N > 1 and any initial population Xq € Cn, 
lim„_>+oo ?T'^^'^(A'^) = 0, both with probability one and in L^(F). 

Proof. Using the notations of Proposition [H let F/v := ■m-jy'' ~ it^^n • From As- 
sumptions (Al) and (A2), one deduces that E{Fiy) < +oo. Then, by Proposi- 
tion [U one has that, for all n > un, ^{n~^d{X^)) < E{F]^)/n, so that conver- 
gence to in L^(P) is proved. Moreover, for any t > 0, Proposition [T] yields that 
T.n>UN nn'^d{X^,) >i)< En>UN Pi^N > m) < EiFN)/i, so that convergence 
to P— a.s. follows from the Borel-Cantelli lemma, since l can be taken arbitrarily 
small. □ 

3.2. Monotonicity properties. The following lemma states a key monotonicity 
property of our branching-selection mechanism. 

Lemma 1. For all 1 < Ni < N2, and fii € CjVi, ^2 S such that -< H2, there 
exists a pair of random variables {Z^, Z'^) taking values in Cn^ x Catj, such that: 
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• the distribution of Z"^ fori = 1,2 is that of the population of particles obtained 
by performing one branching- selection step (with Ni particles) starting from 
the population fii; 

• with probability one, ^ Z"^. 

Proof. Consider an i.i.d. family (^j,j)jG[i,iV2], j=i,2 with common distribution p. For 
k = 1,2, write fik = Y.i^=i^x^{k), with xi{k) > ... > XN^i^)- Then let Tk := 
=1,2 ^Xi{k)+ei,p and define Z^ as being formed by the A''^ rightmost particles 
in Tfc. From the assumption that fii ~< fi2, we deduce that Xi{l) < Xj(2) for all 
i £ Jl, A'^il, whence the fact that Xi{l) + Sij < Xi{2) + Eij, for all 1 < i < A^i 
and j = 1,2. It is easy to deduce that Ti -< T2, whence Z^ -< Z^. The conclusion 
follows. □ 

An immediate corollary is the following. 

Corollary 2. For all 1 < Ni < N2, and /ii G Cjv-^, fi2 S Catj such that ni -< /i2, 
there exists a coupling {Z\, Z^)n>o between two versions of the branching-selection 
particle system, with Ni and N2 particles respectively, such that Zq := fii, Zq := ^2, 
and Z\ Z^ for all n > 0. 

Proposition 2. There exists VMip) G ^ such that, with probability one, and in 
L^(F), for any initial population Xq E Cat, 

lim n~^minX^= lim n~^maxX;^ =vn{p)- 

Proof. Note that, in view of Corollary [H if either of the two limits in the above 
statement exists, then the other must exist too and have the same value. Moreover, 
owing to the translation invariance of our particle system (the dynamics is invariant 
with respect to shifting all the particles by a translation on M) , and to Corollary [U 
we see that it is enough to prove the result when Xq = N5q. The idea of the 
proof is to invoke Kingman's subadditive ergodic theorem (see e.g. [12]), using the 
monotonicity property described by Lemma [TJ 

Consider an i.i.d. family (e£,ij)£>o,iG[i,Ar], j=i,2 with common distribution p (the 
index £ will be used to shift the origin of time when applying Kingman's theo- 
rem). For all £ > 0, denote by (W/^^) fc>o the branching-selection system starting at 

:= N5q and governed by the following steps. For /c > 0, write W^^ = Ei^i ^^i^ 
with xi > . . . > xjq. The population Ti^k derived from W/^ by branching is then 

defined by T^^k ■= E(i,j)eIi,JVlx{i,2} Then, 14^^^^^^ is obtained from T^^k 
by keeping only the A'^ rightmost particles. 

Observe that {W^^) „>o has the same distribution as {X^)n>o. Moreover, the 
argument used in the proof of Lemma [1] shows that 

(5) for all n,m > 0, maxVFj^,,.^ < maxW,]^ -|- maxVF^^. 

Indeed, it is enough to note that ([5]) compares the maximum of two populations 
obtained by performing m branching-selection steps coupled as in LemmalU starting 
respectively from Wq^^ (for the l.h.s.) and from NS^^y^N (for the r.h.s.). 
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Moreover, it is easily seen from the definition that, for each d > 1, the random 
variables {W^^)n>o form an i.i.d. family, and that the distribution of (W^/^,)fc>o 
clearly does not depend on i. One can then check from the definition that the 
following inequality holds |maxWo^„| < Efc=d E(i,j)G[i,Aflx{i,2} Using as- 

sumptions (Al) and (A2), it is then quite clear that there exists ip > —oo such that, 
for all n>0,ipn< E{W^J < +00. 

We conclude that the hypotheses of Kingman's subadditive ergodic theorem hold 
(see e.g. [H]), and deduce that lim„_^+oo uiaxX^ exists both a.s. and in L^(P), 
and is constant. □ 

Proposition 3. The sequence {vN{p))]\f>i is non- decreasing. 

Proof. Consequence of the fact that, when A^i < -^i^o -< ^260, and of the 
monotonic coupling property given in Corollary [2j □ 

We can deduce from the above proposition that there exists Vooip) such that 
limjy^+ooVNip) = V(x,{p)- A consequence of the proof of Theorem [T] below is that 
Vooip) is in fact equal to the number v{p) := A'(t*), which is finite from our assump- 
tions on p. 

3.3. Coupling with a family of N branching random walks. Let 

(BRWj)jg|i TV]) denote independent copies of a branching random walk BRW 
as defined in Section [2l Each BRWj thus consists of a binary tree Tj and a map 
^i. For 1 < i < A^, and n > 0, remember that Ti(n) denotes the set of vertices 
of Tj located at depth n, and define the disjoint union := Ti(n) U • • • U TAr(n). 
For every n, fix an a priori (i.e. depending only on the tree structure, not on the 
random walk values) total order on ■ We now define by induction a sequence 
(G^)„>o such that, for each n > 0, is a random subset of with exactly N 
elements. First, set Gq := Tg'^. Then, given n > and G^, let denote the 
subset of Tj^i formed by the children (each with respect to the tree structure it 
belongs to) of the vertices in G^. Then, define G^_^_l as the subset of formed by 
the A^ vertices that are associated with the largest values of the underlying random 
walks (breaking ties by using the a priori order on T^). Now let denote 
the (random) empirical distribution describing the values taken by the $jS on the 
(random) set of vertices G^ . The sequence (X^)n>o has the same distribution as 
{X^)n>o, when started from := N6o. Thus, we can take for our reference prob- 
ability space (O, P) the one on which BRWi, . . . , BRW at are defined, and let X^ 
be equal to the empirical distribution associated with the subset G^ , and so obtain 
a coupling between {X:^ )n>o (with X^ = N60) and the A^ branching random walks 
BRWi,...,BRWjv. 

4. Results on the branching random walk killed below a linear 

space-time barrier 

Let us start with the following definition, adapted from [16]. Given i; G M and 
m > 1, we say that a vertex u € BRW is (m,w)— good if there exists a finite 
descending path u =: uo,ui, . . . , Um such that ^{ui) — ^{uq) > vi for all i G [[0, mj. 
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Similarly, we say that u is (oo,f)— good if there exists an infinite descending path 
u =: uo,ui, . . . such that ^{ui) — ^{uq) > vi for all i G |0,+oo|I. 

With this terminology, the main result in [T3j can be stated as follows, remem- 
bering that v{p) = A'{t*) and x{p) = ^t*A."{t*). 

Theorem 2. (Theorem 1.2 in [I3]J Let p{oo,e) denote the probability that the root 
of BRW is {oo,v{p) — e)—good. Then, as e goes to zero, 



We shall need a result which, although not stated explicitly in [13], appears there 
as an intermediate step in a proof. 

Theorem 3. (Proof of the upper bound part of Theorem 1.2 in [13jJ Let p(m,e) 
denote the probability that the root of BRW is {m,v{p) — e)—good. For any < (3 < 
x{p), there exists 9 > such that, for all large m. 



One should also consult the papers [T7] for an approach of these results based 
on (mathematically non-rigorous) theoretical physics arguments. See also the dis- 
cussion in Section [H 



The arguments used here in the proof of the lower bound, combine ideas from the 
paper |16] by Pemantle, which deals with the closely related question of obtaining 
complexity bounds for algorithms that seek near optimal paths in branching random 
walks, and the estimate on p{m,e) from the paper [13], by Gantert, Hu and Shi. In 
fact, the proof given below is basically a rewriting of the proof of the lower complexity 
bound in [16j in the special case of algorithms that do not jump, with the following 
slight differences: we are dealing with N independent branching random walks being 
explored in parallel, rather than with a single branching random walk; we consider 
possibly unbounded random walk steps; we use the estimate in [13] instead of the 
cruder one derived in |16j . 

We start with an elementary result adapted from |16| . 

Lemma 2. (Adapted from Lemma 5.2 in [16J.J Let vi,V2 € M 6e such that vi < V2, 
n > 1, m € [[1,?T-]], K > 0, and let =: xq, . . . ,Xn be a sequence of real numbers 
such that Xi-^-l — Xi < K for all i € [0, n — Ij. Let I := {i £ |0, n — m|; Xj — Xi > 



viU - i) for alljGli,i + mj}. If Xn > V2n, then #1 > -K/{K-vi). 



Since Lemma [2] admits so short a proof, we give it below for the sake of complete- 
ness, even though it is quite similar to that in jl6]. 




Proof of Lemma{^ (Adapted from [16j.) Consider a sequence =: xq, . . . , as in 
the statement of the lemma. Let then tq := 0, and, given < n, define inductively 





5. The lower bound 
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Tj_|_i := inf{j € |rj + 1, n|; xj < x^-, + f i(j — Tj) or j = Tj + m}, with the convention 
that inf = n + 1. Now "color" the integers A; G |0, n — 1], according to the following 
rules: if x^^-^ > Xr^ + t;i(rj+i — Tj) and Tj+i < n, then Ti, . . . , Tj+i — 1 are colored 
red. Note that this yields a segment of m consecutive red terms, and that Tj then 
belongs to /. Then color in blue the remaining integers in [[0,n — Ij. 

Let Vred (resp. Vbiue) denote the number of red (resp. blue) terms in [0, n — Ij. 
Then decompose the value of x„ into the contributions of the steps x^+i — x^ such 
that k is red, and such that k is blue, respectively. On the one hand, the contribution 
of red terms is < KVred- On the other hand, the contribution of blue terms is 
< yuue ^ vi + Km, where the m is added to take into account a possible last 
segment colored in blue only because it has reached the index n. Writing that 

n = yred + yUuei We deduce that V2n < KVred + Vi{n — Vred) + Km, so that V^ed > 

n — Km/{K — vi). Then use the fact that at least Vred/fn terms belong to 
/. ' □ 

In [16] , the result corresponding to our Lemma O and an estimate of the type 
given by Theorem [3l are used in combination with an elaborate second moment 
argument. In the present context, the following first moment argument turns out to 
be sufficient. 

Proof of the lower hound part of Theorem\^ Assume that Xq = N6q. Let /3 > 
and > be as in Theorem [3j Then let A > 0, and define 



m 



03/2 



;i + A)(log7V)^^" 



(X(p)-/3)V2 

and e := O/m"^/^, so that, by Theorem [31 

(6) p{m, e) < Ar-(i+^) for ah large N. 

Then let < 7 < 1 and define V2 := v{p) — (1 — 7)e and vi := v{p) — e. 

Let also n := [N^] for some < ^ < A. Now consider k > 0, and let K := 
Klog{2Nn). Consider the maximum of the random walk steps performed during the 
branching steps of {X^ )k>o between time and time n. There are 2Nn such steps, so 
that, by assumption (A2), there exists a value of k such that the probability that this 
maximum is larger than or equal to K is less than {2Nn)~'^^^^ for all large enough N. 
Now denote by Bn the number of vertices in Gq U • • -UG^ (see Section [33]) that are 
(m, v{p)—e)—good (each with respect to the BRW, it belongs to). Observe that, with 
our definitions, for N large enough, ^rfj"^ ~ K/{K — vi) > 0. As a consequence, 
using Lemma [51 we see that, for N large enough, the event maxX^ > V2n implies 
that either there exists a random walk step between time and n which is > K, or 
Bn > 1. Using the union bound and the above estimate, we deduce that 

(7) P (maxX^ > V2n) < (2A^n)"2°°® + P(5„ > 1). 
On the other hand, -B„ can be written as 

(8) Bn:= Yl H u is {m,v{p) -e)-good)l{u eG^ U---UG^). 

ueiiU-'-ulN 
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Now observe that, by definition, for a vertex u at depth £, the event u € Gq U 
• • • U is measurable with respect to the random walk increments performed 
at depth at most i, that is, the family of random variables ^i{w) — ^i{v), where 
i G p,A/']], v,w £ Tj, w; is a child of v (with respect to the tree structure of Tj), 
and w,v are both located at a depth < ^ in Tj. On the other hand, the event that 
u is {m,v{p) — e)— good is measurable with respect to the random walk increments 
performed at depth at least £, that is, the family of random variables ^i{w) — ^i{v), 
where i G [1, A^|, v,w £ Tj, w is a child of v, and w, v are both located at a depth 
>im 1i. 

As a consequence, the two events {u G Gq U- • -UG^ } and {u is (m, v{p) — e)— good} 
are independent. Since the total number of vertices in Gq U • • • U G^ is equal to 
N{n + 1), we deduce from ([SD and ([8]) that E(B„) < N{n + l)N-^^+^\ Using 
Markov's inequality, we finally deduce from ([7]) that 

(9) P(maxX^ > V2n) < (2A^n)-2°°® + (n + 1)N-^. 

Now start with the obvious inequality, valid for all t, exp(tmaxX^) < 
«eXi(n) exp(t$j(u)). Taking expectations, we deduce that E(exp(t max 
A^2" exp(nA(t)). Using the definition of t* and v(p), we then obtain that 

(10) E(exp(t*(maxA:^ - v{p)n)) < N. 
Using (jlOp . we deduce thalQ, for all 6 > 0, and all large enough n, 

(11) E [maxX^l(max A:^ > {v{p) + b)n)] < A^exp {-^t*bn) (1 + \v{p)\n). 
Now observe that, by definition, E(n~^ maxX,^) is bounded above by 

V2 + (vip) + 6)P(maxX^ > V2n) + n^^E [maxX^l(maxX^ > {v{p) + b)n)] . 

Choosing a 6 > , we deduce from ([9]), pT]) . and the definition of V2, that, for all 
large enough N, 

E(n-VaxX^) < {v{p) - {1 - -/)e) + o{{log N)'^). 

Using subadditivity (see the proof of Proposition [2]), we have that VNip) < 
E(n~^ max AT^), and we easily deduce that 

vn{p) < {v{p)-{l-j)e) + o{{logN)-^). 

Now remember that, as N goes to infinity, e ~ ^ff^^^f (log N)^'^. Since the above 
estimates are true for arbitrarily small /3, A and 7, the conclusion follows. □ 



Here are the details. Let M := maxX^ — v{p)n. From the fact that, for all large enough x, 
X < exp(t*a;/2008), we deduce that E(M1(M > bn)) < Eexp{t*M - ^t*bn) for aU large enough 
n. Similarly, E(i;(p)nl(M > bn)) < \v{p)\n¥, exp{t* M — t*bn). The result follows from summing 
the two inequalities above and applying (fTo)). 
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6. The upper bound 

The proof of the upper bound on Voo{p) — vn{p) given in [4j was in some sense a 
rigorous version of the heuristic argument of Brunet and Derrida according to which 
we should compare the behavior of the particle system with particles, to a version 
of the infinite population limit dynamics suitably modified by a cutoff. The proof 
given here relies upon a direct comparison with branching random walks, using the 
fact that, above the threshold induced by the selection steps, the behavior of our 
branching-selection particle system is exactly that of a branching random walk. 

Consider < A < 1 and let e := (^(^i_x^og j^^^ ■ With this choice of e, as goes to 
infinity. Theorem [2] yields that 

(12) p(oo,e) = iV-(i-^)+°W. 

Let us now quote the following result, which is a consequence of Theorem 2 Section 
6 Chapter 1 in [1]. 

Lemma 3. Let (M„)„>o denote the population size of a supercritical Galton-Watson 
process with square-integrahle offspring distribution started with Mq = 1. Then there 
exist r > and (j) > 1 such that, for all n > 0, 

P{Mn > (A") > r. 

Let R be such that R < v{p) and p{[R, +oo)) > 2/3. Consider a Galton-Watson 
tree whose offspring distribution is a binomial with parameters 2 and p{[R, +oo)). 
The average number of offspring is thus equal to 2p{[R,+oo)) > 4/3 > 1 with our 
assumptions. In the sequel, we use the notations r and (j) to denote the numbers 
given by Lemma [3] when we use this offspring distribution. 

Now, let Sn := f^^^] + 1, consider < < 1, and define m := |"MPlzj5)£iv "| ^^^^ 
n := m + S]\f. Let u denote a vertex at depth m in a branching random walk BRW, 
and assume that ^{u) > {v{p)—e)m. Consider the probability that, conditional upon 
the values of $ on the vertices located at depth at most m, there are at least cf)^^ 
distinct descending paths u =: Um, ■ ■ ■ ,Un starting at u and satisfying Uj+i — Ui>R 
for all i € |m, n — Ij. Lemma [3] above shows that this probability is > r. Moreover, 
with our definition of m and n, and our assumption on the value of ^{u), any such 
descending path has the property that $(ni) > {v{p) — e(l -|- r]))i for all i G |Im,n|. 
We conclude that the probability that there exist at least (j)^'^ distinct descending 
paths of the form root = uq, . . . ,Un such that ^{ui) > {v{p) — e(l -|- rj))i for all 
i € |IO,n|, is > p{m,e)r. 

Now define A as the event that, for all j € Jl,A^|, BRWj does not contain 
more than (f)^^ distinct descending paths of the form root = uq,. . . ,Un such that 
^j{ui) > iv{p) - e(l + r]))i for all i G [O, nj. Using the fact that BRWi, . . . , BRW^v 
are independent and the above discussion, we see that 

IP(^) < [I- p{m,e)rf . 

Using (fT2]) . the obvious inequality p{m,e) > p{oo,e), and the fact that 1 — x < 
ex.p(—x) for all x, we deduce that 

(13) P(^) <exp(-iV^+°W). 
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Let S := e(l + rj). Define the event B := {min(X^) < {v{p) — 6)k for all k € 
p,n|}, and assume that B D occurs. From the definition of the selection mech- 
anism, we conclude that there must be at least (/>*^ distinct vertices in the set G^, 
which is a contradiction since > N. As a consequence, B (1 A'^ = 9, so that 
B cA. 

From (1131). we thus obtain that 



(14) P(S) < exp(-Af^+°(i)). 
To exploit this bound, we use the following result. 

Proposition 4. With the previous notations, for all N large enough, 

vn{v) > {v{p) -S)- \v{p) - 6\nF{B) - nE(IG„|l(S)), 

where 0„ is the minimum of 2nN i.i.d. random variables with distribution p. 

Proof. We re- use the coupling construction given in the proof of Proposition [21 and 
assume that {X^ )n>o is defined using this construction by the identity := W^^. 
Start with Fq := and Jq := 0, and i := 0. Given i > 0, Fj and Jj, let Lj+i := 
inf{fc € min(VF]^^) > {v{p) — 6)k}, with the convention that inf := n. Then 

let Fj+i := Ti + Lj+i, and let Jj+i := Jj + min(Wj^^^^^^). 

Using an argument similar to the proof of Lemma [H it is then quite easy to deduce 
that, a.s., 

(15) for all i > 0, minWo^r, > Ji- 

Observe that the sequence (Fj+i — Fi)j>o is i.i.d., and that the common distri- 
bution of the Fj_|_i — Fj is that of the random variable L defined by L := infj/c G 
[[l,n|; m.m{X^) > {v{p) — 5)k^, with the convention that inf := n. Similarly, the 
sequence (Jj+i — Ji)i>o is i.i.d., the common distribution of the J^+i — Jj being that 
of minX^. 

From the law of large numbers and Proposition [21 we have that, a.s., 
limj_^_^oo min = vn{p)K{L), while the law of large numbers and ^S) imply 

that liminfj_j.+oo minX]?^ > E(minX^). We conclude that V]\j{p) > ^^"e(L)'^ ^ • 
Now, let G„ denote the minimum of all the random walk steps performed by the 
branching-selection system between time and n. 

By definition we have that minX^ > {v{p) — 6)L1{B'^) + L©nl{B), so that 
E(minX^) > (vip) - (5)(E(L) - E(L1(5))) + E(Le„l(S)). Using the fact that 1 < 

L<n, we obtain that ^^'^'("^^'^^^ > {v{p) - 6) - \v{p) - 6\nF{B) - nE{\en\l{B)). □ 

Proof of the upper bound part in Theorem Ul In view of Proposition [H we deduce 
that vn{p) > {vip) - (1 + ??)e)(l - nF{B)) - nE(|e„|l(5)). Bounding above |e„| 
by the sum of the absolute values of the 2nN corresponding i.i.d. variables, 
and using Schwarz's inequality thanks to Assumptions (Al) and (A2), we deduce 
that E(|0„|l(i?)) < 2nNCF{By/'^ for some constant C (depending only on p). 
From (fill) and the definition of n, we deduce that, as N goes to infinity, nF{B) and 
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nE(|6„|l(B)) are o((logiV) so we obtain that 

VNip) > v{p) - ^^M(i_t^(iogAr)-2 + o((logiV)-2). 

(i-xy 

Since A and ij can be taken arbitrarily small in the argument leading to the above 
identity, the conclusion follows. 

□ 

7. The Bernoulli case when 1/2 < a < 1 

In the Bernoulli case p = a6i + (1 — a)6o, with 1/2 < a < 1, Assumption 
(A3) breaks down, and the behavior of the particle system turns out to be quite 
different from Brunet-Derrida, as stated in the following theorems. Note that, when 
1/2 < a < 1, v^{p) = 1. 

Theorem 4. For a = 1/2, there exists < c*(p) < c*{p) < +oo such that, for all 
large N, 

(16) c^{p)N-^ < 1 - vn{p) < c*{p)N-\ 

Theorem 5. For a > 1/2, there exists < d*{p) < d^{p) < +oo such that, for all 
large N, 

(17) exp(-4(p)iV) < 1 - vn{p) < exp{-cl*{p)N). 

7.1. Lower bound when a = 1/2. It is easily checked that, for all m > 0, 
the number of particles in the branching-selection system that are located at po- 
sition m after m steps, that is, X^{m), is stochastically dominated by the total 
population at the m— th generation of a family of independent Galton- Watson 
trees, with offspring distribution binomial(2, 1/2). This corresponds to the crit- 
ical case of Galton- Watson trees, and the probability that such a tree survives 
up to the m— th generation is < cm~^ for some constant c > and all large 
m. As a consequence, the union bound over the N Galton- Watson trees yields 
that, for large enough m, ¥{X^{m) > 1) < cNm~^. On the other hand, we 
have by definition that Emax(X^) < mW{X^{m) > 1) + (m - l)P(X^(m) = 0). 
Choosing m := AN, where A > 1 is an integer, we deduce that, for large A^, 
m~^Emax(X^) < 1 — -xn{1 — c/A). Using subadditivity (see the proof of Proposi- 
tion[2]), we have that V]\i{p) < K{m^^ max.X^). The lower bound in (I16p follows by 
choosing A > c. 

7.2. Upper bound when a = 1/2. Given m > 1, define U := inf{n G 
[[l,m]]; X^ (n) < 2N/3}, with the convention that inf0 := m. Observe that 
minX^ > U — 1, since, by definition, X^_^{U — 1) > 2N/3, so that, after the 
branching step applied to X^_^, the number of particles whose positions are >U — 1 
must be > 2 X 2A^/3, whence > A^. 

Using an argument similar to the proof of Proposition HI we deduce that 

, , 1 



14 



JEAN BERARD, JEAN-BAPTISTE GOUERE 



The lower bound in ()16p is then a direct consequence of the following claim. 
Claim: for small enough e > 0, with m := [eiVj, there exists c(e) > such that 
E(f/) > c(€)N for all large A^. To prove the claim, introduce for every x G N the 
Markov chain {V^)k>o defined by the initial condition Vq := x, and the following 
transitions: given Vq , . . . ,VI^ , the next term V^_^_i is the minimum of A'^ and of 
a random variable with a binomial(2V^^, 1/2) distribution. Clearly, the sequences 
{Vj^)k>o and {X^ {k))k>o have the same distribution. Moreover, given two starting 
points x,y gN such that x <y, one can easily couple {V^)k>o and (V^^)fc>o in such 
a way that < for all k > 0. As a consequence, choosing xn ■= [3-^/4], 
we see that U stochastically dominates the random variable T defined by T := 
inf {n E V^^ < 2N/3} (again with inf := m), so that F{U = m) > P{T = 

m). 

Now let us define yet another Markov chain {Zk)k>o by Zq := xjy and the 
following transitions: given Zq, . . . , Z^, the next term Z^+i is a random variable 
with a binomial(2Zfc, 1/2) distribution. Clearly we can couple {Zk)k and (V^^^)fc 
so that they coincide up to one unit of time before the first hitting of [A'", +cx)[[. 
As a consequence, the two events Ai := {sup^g|o„^j — [3A^/4J| < A^/16} and 
A2 := {sup^g|o_m] — [3-^/4] I < A^/16} have the same probability. Observing 
that {Zk)k>o is a martingale, we can use Doob's maximal inequality to prove that 
P{A1) = P{A^2) < E{Zm - [3N/4\)'^{N/I6y^. Then, it is easily checked from 
the definition that E{Z'^^^\Zk) = Z| + Zk/2 for all A; > 0, and, using again the 
fact that {Zk)k>o is a martingale, we deduce that E{Zm — [3A^/4J)^ < mN/2. 
As a consequence, we see that, choosing e > small enough, we can ensure that 
-P(^i) < 1/2008 for all large A^. Since, by definition, Ai implies that T = m, 
we finally deduce that, for such an e, and all A^ large enough, we have that 
P([/ = m)> P{T = m)> 2007/2008. The conclusion follows. 

7.3. Upper and lower bound when 1/2 < a < 1. As for the lower bound, observe 
that the probability that all the 2A^ particles generated during a branching step 
remain at the position from which they originated is (1 — a)2^, so that E(maxA'^) < 
n(l - (1 - af^). As for the upper bound, observe that, starting from A^ particles 
at a site, the number of particles generated from these during a branching step 
and that perform +1 random walk steps has a binomial(2A'^, a) distribution, whose 
expectation is 2aN, with 2a > 1. Using a standard large deviations bound for 
binomial random variables, we see that the probability for this number to be less 
than N is < ex.p{—cN) for some c > 0. Using superadditivity E(minX^) (derived 
in exactly the same way as the subadditivity property of E(maxX^), see the proof 
of Proposition [2]) , it is easy to deduce that E(min X^) > n(l - exp(-cAr)). The 
result follows. 

8. Discussion 

This section contains a discussion whose goal is to provide a self-contained qual- 
itative explanation of the (logA^)~^ order of magnitude appearing in Theorem [H 
Most of the discussion consists in explaining the e~^/^ scaling of log p(oo, e), and of 
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log p{m,e) when m oc e^^^^, in a way that is (hopefully) less technically demand- 
ing than the proofs presented in [13], although we follow the proof strategy of [13] 
rather closely. Note that the discussion here deals mostly with the order of magni- 
tude of terms, not with the precise value of the constants as in [13]. For the sake of 
readability, some of the arguments are only discussed in a quite informal way. 

8.1. Asymptotic behavior of p{oo,e) and p{m,e). 

8.1.1. Connection between p{oo, e) and p{m, e). A first remark is that the asymptotic 
behavior of p{oo, e) can be connected with that of quantities of the form p{m, e) under 
appropriate conditions. One obvious inequality, valid for all m > 0, is the following 

(19) p{m,e) > p{oo,e). 
In the reverse direction, we have the following. 

Proposition 5. There exist R < v{p) — 1, > 1, r > and c > 0, depending only 
on p, such that, for all m > 0, and all < e < 1, the condition 

(20) (p'^pim, (1 - a)e) > c 
implies that the following inequality holds 

(21) p{oo,e)>^pim,{l-a)e), 

where a and e are arbitrary numbers satisfying < a < 1 and e > 0, and 

aem 



^v{p) — e — R^ 

The proof of the above proposition uses the following elementary lemma. 

Lemma 4. Consider a Galton- Watson process with offspring distribution Q. If 
there exists a > 1 such that a x Q{[a, +oo[) > 2 log 2, then the survival probability is 
larger than or equal to Q([a, +oo[)/2. 

Proof of Lemma^ Let g{s) := X^^^ Q{k)s^ for s G [0, 1[. By coupling, it is enough 
to prove the result under the additional assumption that only the values and 
a have non-zero probability with respect to Q, so we may assume that g{s) = 
l — Q{a) + s"'Q{a). Since aQ{a) > 1, we have a super-critical Galton- Watson process, 
and, from standard theory, we know that the extinction probability d of the process 
is the unique solution in [0, 1[ of the equation g{d) = d, with g{s) > s for s €]0,(i[ 
and g{s) < s for s €.]d, 1[. Our assumption that a x Q{a) > 2 log 2 easily yields the 
fact that g{l — Q{a) /2) < 1 — Q{a)/2, whence the fact that d must be < 1 — Q{a)/2. 
The result follows. □ 

Proof of Proposition\M Consider the values of R, 4>, r defined in Section [H in the ar- 
gument following Lemma [31 Then consider a descending path root = uq, . . . , Um+q G 
T such that ^{ui) > {v{p) - (1 - a)e)i for ah i G [[0, m], and ^(ui+i) - ^(ui) > R 
for all i € [[m,m + q — 1|. We see from the definition of q that <J>(u,) > {v{p) — e)i 
for all z G do, m + g|. 

We now define a Galton- Watson branching process of vertices of T in which, for 
all n, the n— th generation of the process is formed by vertices in %{{m + q)n). First, 
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the zero-th generation of the process is formed by the root of T. Then, given a vertex 
X G %{{m+q)n) belonging to the n— th generation of the process, the offspring of this 
vertex in the branching process is formed by all the endpoints y of descending paths 
X =: uq, . . . , Um+q := y in T such that ^{ui) — ^{uo) > {v{p) — e)i for all i G |0, m+g|. 
Prom the definition of the branching mechanism of BRW, we see that we have defined 
a Galton- Watson branching process. Now, re-doing the argument following LemmaO 
in Section [H we see that the offspring distribution of this branching process gives at 
least 0^ children with probability at least p{m, (1 — a)e)r. On the other hand, the 
definition of our branching process shows that if it never goes extinct, the root of 
T is {oo,v{p) — e)— good. As a consequence, the survival probability of our process 
is a lower bound for p{oo,e). The result then follows from Lemma [H choosing 
c>(21og2)/r. □ 

We conclude this section by the following remark: using the above results, it is 
possible to deduce the conclusion of Theorem [3] from the conclusion of Theorem [2l 
In [13] , Theorem [3] is in fact an intermediate step in the proof of Theorem [2l so 
our remark does not lead to an alternative way of proving Theorem [3] from first 
principles. However, its interest is to show that, as soon as e^^"^ log p{oo, e) converges 
to some limit, this limit can be approached arbitrarily closely by expressions of the 
form e^/^ logp(m, e), with e = Qjrr?!^ for some large enough constant Q. 

Proof of Theorem\^ from the conclusion of Theorem\M Let 0<a<l,6>0 and 
m > 1. Set e := Q jrr?!'^ . Let i?, (/>, r, c, q be defined as in the statement of Proposition 
[5l As m goes to infinity, we have that log^'^ ~ log(^)(i'(p) — and, by 

Theorem [21 

logp(oo, (1 - Q)e) ~ -x(p)^/^^"^/^(l - aY^I'^m^l'^. 

Therefore, provided that d has been chosen large enough, Inequality (j20p holds for 
large enough m. Given such d and m. Proposition [5] and Theorem [2] yield that 

p{m, (1 - a)e) < (2/r)p(oo, e) = (2/r) exp ' ^^^^^ ^ ""^"^^ ^ ' 

Setting (1 — a)0m~^/^, one obtains that 

p{m,'e) < exp 



x(p)(l-a)+o(l) 



1/2N 



Since a can be chosen arbitrarily small, the conclusion follows. □ 

8.1.2. Strategy and results. Given the results of the previous section, the strategy 
consists in studying the order of magnitude of log p{m, e), when m has the scaling 
form m oc e~" for some n > 0, seeking a value of u such that (|20|) is satisfied. We 
shall see that logp(m,e) oc e~^^^\ with h{u) = u/3 for all < u < 3/2, while the 
integer q in Proposition [5] satisfies q oc e~^"~^^ for all u > 1. For u* := 3/2, the 
identity h{u*) = u* — 1 holds, so that q and log p{m, (1 — a)e) have the same order 
of magnitude and (f20]l can be satisfied. Then (fT9]l and (f2T]l imply that log p{oo, e) oc 
—g-'^C" ) = _£-i/2^ g^j^d (using the fact that h is non-decreasing), h{u) = h{u*) for 
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all u > u* . Thus, the e~^/^ scaling exponent is "explained" by 3/2 being the solution 
of the equation h{u*) = u* — 1. 

We deduce from these results the existence of two distinct regimes for p{m, e) with 
m cx e~": 

• when < n < 3/2, log(p(m, e)) cx —m^^'^; 

• when u > 3/2, log(p(m, e)) cx log(/9(cx), e)) cx —e^^^"^. 

8.1.3. Asymptotics of log p{m,e) with m cx e~", < u < 3/2. Let us now explain 
how to compute h{u), and assume throughout this section that m cx e~", whence 
em cx e^~". A key idea is to perform a change of measure, replacing the step 
distribution p of the BRW by the distribution p defined by (see Section II. 2p 

dp, , exp 



dp exp(A(t*)) 

The mean value of a step with respect to p is now equal to v{p), and, if {Sk)k>o 
denotes a random walk started at Sq := 0, with i.i.d. increments whose common 
distribution is p with respect to a probability measure P, and p with respect to a 
probability measure P, the following identity holds for all k: 

(22) 2''P [{So, ...,Sk)e-]=E [e-**(^'=~''(^')'=)l(So, . . . , 5^) e 

where E denotes expectation with respect to P. The raison d'etre of Assumption 
(A3) is to allow for such a change of measure. 

Remember that p{m, e) is the probability that at least one descending path root =: 
uo,ui, . . . , Um exists in T such that 

(23) for ah i £ [I0,m|, <^{ui) > {v{p) - e)i. 
Observe that, for such a path, either 

(24) for all i G [0, m|, <i>{ui) < v{p)i + e~"/^ 
or 

(25) there exists i G I0,m]] such that <I>(tii) > v{p)i + e""/^. 

Denoting by the (random) number of descending paths satisfying (j23l) and (I24p . 



and by the number of those satisfying (j23]) and (|25]) . we see that 

p(m, e) = Q({H„ > 1} U {A^ > 1}), 

so that, obviously, 

(26) Q(H„ > 1) < p(m, e) < E(H„) + Q(A„ > 1). 
By definition, 

E(H^) = 2"^P( for all i G [0,m], ?;(p)i - ei < Si < v{p)i + e'^/^), 
which rewrites, using (I22p . as 

(27) ^ L-**(5m~»;(p)m)-^^ g^j^ • ^ ^^^^ • _ < 5. < + g-^/S) 
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Since the only paths that contribute to the above expectation have 

v{p)m — em < Sm < v{p)m + e""''^, 

we see that 

(28) E(H„) < e^'^'^P [ for all i G |0, mj, v{p)i - ei < < v{p)i + e-"/^ ] . 
and 

(29) E{Em) > e-**'""^'P [ for all i G JO, mj, v{p)i - ei < S^ < v{p)i + e^^/^ ] , 

Now observe that, under P, (S^— v(p)A;)fc>o is a random walk with centered square- 
integrable increments. Moreover, (e"'"/^)^ << m since u > 0, and, e~"/^ >> em as 
soon as u < 3/2. As a consequence, the usual Brownian scaling for random walks 
yields that when < n < 3/2, 

(30) logP [ yi e [[0,m], v{p)i -ei<Si< v{p)i + e-"/^ 1 ^ _ ^ " = -e-"/^. 

Since e~"/'^ >> em when u < 3/2, we deduce from (p8|) and (p9]) that 

(31) logE(H„) oc -e""/^ for 1< < 3/2. 
As for Am, the union bound yields that 



> ^ _ 

1=0 

whence 



1) < X] Q(3x G T(i); $(x) > v{p)i + e-"/^), 



W > 1) < E > ^b)^ + e""^')- 

For all i G [0, m|, the change of measure shows that 2*P(5.j > v{p)i + e~"/^) 
P [e~**(-5-^{p)»)l(5. > y{^p)i + e-V3)] < exp(-t*e-"/3)^ gg^gj^y deduce that 

(32) - logQ(Am > 1) is at least oc -e""/^. 



We conclude from ([M]), (IMl) and ((26]) that 

(33) - logp(m, e) is at least oc e'^/^ for < n < 3/2. 

Using the same kind of argument that led to (|3T]) . but working a little more (we 
omit the details), it is possible to show that 

(34) logE(H^) is at most oc e""/^ for < u < 3/2. 
Then, we can use the classical second moment inequality: 

> 0) , 

to deduce from ([26]), ([31]) and ([HD that 

(35) - log p(m, e) is at most oc -e~"/^ for < u < 3/2. 
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As a consequence, we obtain that 

(36) log p{m, e) cx -e""/^ for < tt < 3/2. 

8.1.4. Asymptotics of p{m,e) with m oc e"", u > 3/2, and p{oo,e). The above 
discussion dealt only with rough order of magnitudes (denoted by the oc symbol), 
but a more precise analysis is needed to study the competition between positive and 
negative terms of similar orders of magnitude when u = 3/2. 
For A > 0, let m := [Ae^^/^J and 

/^(A) := limsup — e^/^ logp(m, e), /~(A) := liminf — e^''^ log/9(m, e). 

Prom the monotonicity property: p{m, e) > p{m', e) when m' > m, we deduce that 
A I—)- /^(A) and A i— )• /~(A) are non-decreasing. Similarly, from the monotonicity 
property p(m, e) < p{m, e') when e < e', we deduce that A ^ \~'^/^f+{X) and 
A !->■ A~^/^/~(A) are non-increasing. 

In particular, the fact that there exists some A for which /^(A) is finite (resp. 
positive) implies that /'^{X) is finite (resp. positive) for all A > 0. The same property 
holds for /~. 

Now, rework the bounds in the previous section, replacing the e~"/'^(= e"^/^ since 
u = 3/2) terms in the definition of and A^, by Ae~^/^ (more precision than only 
the order of magnitude of terms is needed in order to deal with the case u = 3/2). 
Consider the analog of the bound (I28p in the present context, and observe that, for 
small e, em ~ Ae~^/^. The Brownian scaling bound then yields the existence of a 
constant c > such that, for small e, 

logP [ Vi G 10, ml v{p)i -ei<Si< v{p)i + Ae-i/2 ] < _c_^^ ^ -cX-^-^''^. 

For small enough A, this term dominates the t*em ~ t*Ae^^/^ term in the exponen- 
tial, so that /~(A) > 0. We deduce that /~(A) > for all values of A > 0. On the 
other hand, it is straightforward to adapt the estimates in the previous section to 
show that /^(A) < -|-oo for all A > 0. We can thus conclude that, when m oc e~^/^, 

(37) logp(m, e) oc — e^"*^/^. 

Now, by Proposition m the asymptotic scaling logp(oo,e) oc — e^"*^/^ is a con- 
sequence of (j37|) . provided that (pOj) is satisfied for u = 3/2 and, at least, large 

enough A. Observe that, on the one hand, q ~ v{p)^R ~ 'v{p)-R • 
other hand, \ogp{m,{l - a)e) > -/+(A(1 - af/'^){{l - a)e)-i/2_ xhe fact that 
A I— > A~^/^/"'"(A) is non-increasing and thus bounded above for large A, implies that 
(p'^p{m, (1 — a)e) » 1 for large enough A, so that (pOj) is indeed satisfied. 

8.2. Deducing the Brunet-Derrida behavior. Broadly speaking, our proof of 
the Brunet-Derrida behavior of branching-selection systems is based on the fact that 
there is a loose equivalence between the following two properties: 

(38) BRWi, . . . , BRWat do not survive killing below a line of slope v — e, 
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and 

(39) vjy < V — e. 

If one accepts this premise, it is then natural to expect the actual velocity shift 
ei^ := v{p) — vn to satisfy 

(40) p(oo,e7v) ocl/iV. 

Indeed, since BRWi, . . . ,BRWAr are independent, p{oo,ej\f) >> ^/N would imply 
that, with probability close to one, at least one of the BRWjS survives killing, while 
p{oo,€]\r) « 1/-^ would imply that, with probability close to one, none of the 
BRWjS survives. Using the asymptotics 

logp(oo,e) x(p)^^^e"^^^ 

it is then easily checked that (I40p imposes the precise asymptotic behavior 

eAf ~ x(p)(logiV)"^. 

To give an intuition of why (p8|) and (|39|) should be related, remember the coupling 
between the branching-selection particle system and BRWi , . . . , BRW ^ described in 
Section [3^ If (f38]l holds, then, loosely speaking, v{p) — e is above the sustainable 
growth speed for a branching system with only particles available, so the maxi- 
mum of the branching-selection system should grow at a speed lower than v{p) — e. 
Conversely, if holds, the population in BRWi, • • • ,BRW7v above the line with 
slope v{p) — e quickly exceeds since the existing surviving particles quickly yield 
many surviving descendants. As a consequence, the threshold in the selection steps 
of the branching-selection system has to be above the line with slope v{p) — e. 

A rigorous formulation of the preceding arguments is precisely what we do in 
Sections [5] and El 

It should be noted that the key time scale over which we have to control the 
particle system to prove both the upper and the lower bound is e~^/^, with e satisfy- 
ing (140 p . that is, a time scale of order (log A'^)^. This is the same order of magnitude 
as the one observed for coalescence times of the genealogical process underlying the 
branching-selection particle system (this question is investigated empirically and 
with heuristic arguments in e.g. [5]). Understanding more precisely the role of this 
time-scale for the dynamics of the particle system certainly deserves more investi- 
gation. 
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